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Free electron approximation for electron states in poly(p-phenylene vinylene) and other conjugated 
systems is proposed. It provides simple and clear analytical expressions for energies of electron states 
and for wave functions. The results are in qualitative agreement with experiments. Our model does 
not contain fitting parameters. We present two examples of the developed theory applications: exact 
calculation of electron energy in magnetic field, and scattering of electrons in copolymer. 
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I. INTRODUCTION 

The valence effective Hamiltonian method is used for 
many years as a theoretical basis for band structure 
calculation of conjugated polymers, see Ref. |^ and lit- 
erature cited therein. However, very simple approx- 
imate calculations can-be performed in framework of 
introduced by Schmidta free [Electron model, see more 
references in book of SalemB. This method is use- 
ful for classification of states of one-dimensional chain 
molecules and aromatic molecules. In ijhe latter case 
it is known as Piatt's Perimeter Mfldela or Perimeter 
Free Electron Orbital Theory(PFEO)El However, arylene 
based polymers, like poly(p-phenylene vinylene) (PPV) 
or poly(l,4-naphthalene vinylene) (PNV), are not purely 
one-dimensional, because they contain branching points. 
The generalization of free electron approach to this case 
is the main purpose of the present work. For the first 
time such calculations were undertaken by Ruedenberg 
and ScherrO and here we develop the free electron net- 
work model starting with symmetry arguments. 

We also would like to consider the basic question why 
there exists significant difference in characteristic energy 
scales of polyacetylene and PPV 7r-electron bands siruc- 
ture. The answer was given in the work of Soos et au. In 
the cited work authors addressed the small energy scale 
of the band structure of PPV to the presence of phenyl 
rings. The bands of PPV are formed from the levels of 
isolated phenyl ring (benzene molecule) . The band width 
is determined by the coupling between the rings, which 
can be obtained by the canonical transformation from 
the coupling between adjacent carbon atoms. 

In the present work we calculated the structure of PPV 
TT-bands by means of constructing of 7r-states from free 
electron wave functions. In our model the proposed in 
Ref. 1^ topological gap is the consequence of the time- 



reversion symmetry of the wave function near the triple 
nodes. Triple node is carbon atom bound to three other 
carbons. Therefore, the structure of 7r-bands of PPV 
is determined by the phases, which are accumulated by 
the wave function having passed from one triple node to 
another. Polymer deformations become unimportant in 
this model. However, the proposed here model is suit- 
able only for organic structures build from conjugated 
carbon chains. Other kinds of organic structures cannot 
be described by this model. 

We obtained PPV 7r-electrons wave functions in sim- 
ple analytical form and their symmetry properties are 
explicitly seen. Symmetry of all PPV states is essentially 
different from the symmetry of benzene molecule states, 
because spatial oscillations period of free electron wave 
function changes continuously with energy. One cannot 
say that this or that PPV electron state originate from 
specific benzene molecule state. For example, almost all 
direct light absorption processes between two highest oc- 
cupied bands and two lowest empty bands are dieale al- 
lowed, that agree with observed experimentallyu'cl four 
absorption bands. 

Effect of magnetic field can be taken in to account ex- 
actly in framework of our model. In the same way one 
can explain roagnetic susceptibility anisotropy of ben- 
zene moleculetj. Magnetic field induces ring current in 
benzene molecule and this effect gives main contribution 
in magnetic susceptibility anisotropy. In PPV magnetic 
field also induces currents in phenyl rings, and they give 
main contribution into PPV diamagnetic susceptibility. 
It is different from the susceptibility of equivalent number 
of isolated phenyl rings, because magnetic field induced 
shifts of electron energies in PPV and benzene molecule 
are different. 

In pure PPV one can distinguish occupied valence band 
and empty conductance band. In order to obtain certain 



concentration of free carriers, holes in the valence band 
and electrons in conductance band, one can dope poly- 
mer with impurities or create structural defects. Our 
model can be used for calculation of electron and holes 
transmission probabilities through structural defect, if it 
is build from conjugated carbon chains. In this work 
we considered naphthalene molecule inserted in PPV in- 
stead of phenyl ring in order to see how strong it reflects 
carriers. 

The accurate mathematical formulation of our model 
is done in next section. The PPV and PNV bands struc- 
ture calculated too. Developed method application to the 
calculation of the diamagnetic susceptibility is given in 
Sec. ni. An example of electron scattering on monomer 
of a different kind embedded into polymer, is considered 
in Sec. |^. In the same section we give the explicit form 
of PPV wave functions. The results are summarized in 



Sec. 0. 



II. ELECTRON STATES OF PPV AND PNV. 

Perimeter Free Electron Orbital Theory (PFEO)Q is a 
good tool for description of benzene molecule 7r-electrons 
states. In this approximation the local potential of the 
atoms is ignored, and it's assumed that electrons can 
move freely around the molecule ring. The electron wave 
function for such state is "0 = exjp{ikr), where r is coor- 
dinate along the circular ring, and k is the wave number. 
The periodic boundary condition results in kL = 27rn, 
where L is the perimeter of the ring that is approximately 
8.3A. So this theory gives sequence of benzene 7r-electron 
levels 
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0,1,2,. 



(1) 



Since we have one 7r-electron per each carbon atom, these 
six electrons occupy two lowest energy levels( levels with 
n > 1 are twice degenerate). 

The big surprise is that the results of the above model 
are in good agreement with experiment in spite of ben- 
zene molecule complex structure. For example, carbons 
create periodic potential V{r) along benzene ring perime- 
ter. The first order correction to the energy of n-th state 
due to the presence of this potential is given by its matrix 
element (n|y|n). This matrix element will be different 
from zero if the period of the potential is equal to mul- 
tiplied by integer number half the wave length of n-th 
state wave function. The period of V{r) is i/3 due to 
the alternation of the C-C bonds lengths, and the matrix 
element (n|T/|n) is not zero for the state with n = 3, 6, . . ., 
but other states don't feel this potential. Therefore, the 
application of Eq. (|l|) is justified for low energy states 
n< 2. 

The main idea of this work is to extend the above ap- 
proach to PPV, see Fig. |l|, which can be considered as a 
sequence of benzene molecules connected via short chain 




FIG. 1. Molecular structures of poly(p-phenylene viny- 
lene) and poly(l,4-naphthalene vinylene). Schematic dia- 
grams show choice of coordinates along carbon chains for our 
model. The short vertical lines mark the origins of corre- 
sponding coordinates. 



of two additional carbons. In the framework of the con- 
sidered model, electron runs along this short chain unless 
it reaches the phenyl ring. Here electron has some proba- 
bility to come into the ring or to go back. The same pic- 
ture repeats itself near the exit from the phenyl ring. We 
will obtain the valid state of PPV if electron wave func- 
tion has the same amplitude after passing one monomer. 
Therefore, each unit cell of the polymer is modeled by the 
three bonds, as it is shown in Fig. 0, where the length of 
each bond is L/2. The electron can go from one bond to 
another at the nodes, which are marked by the open and 
filled circles in Fig. ||. The nodes correspond to carbon 
atoms bound to three other carbons. 

Let's imagine the free electron wave is propagating 
along the bond. When the wave reaches the node the 
wave has amplitude r to be reflected and amplitude t to 
go to each remaining bond, see Fig. y. The amplitudes 
to go left or right will be the same if we assume that the 
node has C^y symmetry and the electron is in the ground 
state of the bond confinement potential. More rigorously, 
let's introduce coordinates a;, y, z along three bonds at- 
tached to one node, black circle in Fig. IT], with origins 
at the middle of the corresponding bond. The electron 
wave functions in these bonds will be 



Be 



— ikx 






(2) 



Each wave leaving the node has to be linear combination 
of all the three waves coming in to this node, 

E^-ikLI'^ = rFf^^^l^ + tAf^^^l'^ + toe^^l'^ , (3) 

where phase kL/A is acquired by electron having passed 
from the middle of the bond to the node. 

The amplitudes t and r are complex numbers defined 
by the properties of the node. These amplitudes are not 
completely independent due to time reversion symmetry, 
see Fig. 0. The simple calculation results in two condi- 
tions 



' + 2|t|2 = l. 



(4) 



Therefore, two complex amplitudes t and r are defined 
by two real parameters, which can be chosen as following 




(5) 

Here < VF < 8/9 and < 6* < 2tt. The time rever- 
sion symmetry leads to some unusual property of "triple" 



node transmission and reflection amplitudes. The abso- 
lute values of these amplitudes satisfy the conditions 



< |t| < 2/3 , 1/3 < |r| < 1 



(6) 



rather than the usual < |r| < 1, and < |t| < 1. 
Therefore, symmetric "triple" node always reflects elec- 
trons with the probability larger than 1/9, and this is 
consequence of the time reversion symmetry. 

Our purpose is now to find the states in PPV. Since the 
polymer is periodic structure, the electron wave function 
is characterized by quantum number cj). When electron 
pass one monomer, the wave function acquires this ad- 
ditional phase. Let say that considered previously two 
bonds, with coordinates y and z along them, join a new 
node, blank circle in Fig. |^, and make a loop. If the coor- 
dinate along third bond attached to this node is x' then 
the wave function in this bond will be 



J'l' 



yle'" 



Be 



-ikx' \ 



(7) 



Application of the transmission and reflection rules to 
this node gives us three additional equations 

Ae-'^^l^+^'l' = rBe'''^'^+"^ + tCe'''^'^ + tEe'^^l^ , 

p^-ikL/4 ^ rE^ikL/A ^ ^Q^^kL/i+i4, ^ tCe'^^'^ . (8) 

Six equations, Eqs. (||,||), will be compatible if the ma- 
trix determinant of their coefficients is zero. This deter- 
minant has a following form, as a function of quantum 
number (j) and energy, which is taken into account by 
phase factor exp{ikL), 



^ — ikL 



-{r- tf 



{r + 2tf{r-tf 



— ikL I ^ — 2ikL 



hr^ + 2rt + f + Af cos{(j))\ e-'^^ + e 



(9) 



The equation Q — Q gives the dispersion relation for 
electron states in the polymer. It is a cubic equation for 
exp(z/cL) and its three roots can be represented in the 
following form 

cos(/cL - 61) = (9/4)Ty - 1 , sin(fcL - g) > , 

(10a) 

cos(/cL -&) = W[blA + cos((/))] - 1 , sin(fcL - 6) > Q , 

(10b) 

cos(A:L -9) = W[b/A + cos{cj))] - 1 , sin(fcL -9) <Q , 

(10c) 

where fc > due to our choice of the reflection and trans- 
mission events. The energies of the states are given by 



E = 



2m ' 



(lOd) 



it 



^Z) 



tt 



FIG. 2. The wave coming to the node has the reflection amplitude r, and transmission amplitude to go to each other branch 
i, see left figure. Inversion of time gives two conditions, see right figure. Amplitude of the wave going back has to be one and 
we obtain 1 = tt^ + tt'' + rr^ . Amplitudes of the waves passing to the branches has to be zero and therefore = tr^ + tt^ + rt'' . 



where to is a free electron mass. 



Equation (10a) gives a sequence of isolated levels that 
corresponds to waves trapped in rings. One can easily 
verify that Eqs. (y,||) give for t his case A = B = 0, and 
C = —E, D — —F. Equations ( 10b ,10q) give a sequence 
of bands and envelop wave function is proportional to 
e*'^'^, where v is number of monomer in polymer chain. 
The amplitudes A,B,C = E^D = F can be obtained 
from Eqs. (pH), where k ha s to be cons ider ed as a func- 
tion of <j) given by Eq. ( |lOb| ) or by Eq. (10c). 

The comparison of the energy dispersions with numeric 
results of Ref. |] shows that 9 and W are energy indepen- 
dent and take the limit values: « 0, VF « 8/9. These 
numbers mean that electron does not acquire additional 
phase after passing triple node, and transmission of triple 
node is as good as it is allowed by time reversion sym- 
metry. Therefore, 6 and W take their limit values and 
should_not be considered as fitting parameters. For these 
valuealj transmission and reflections amplitudes are real 
and \t\ takes its maximal value, t = 2/3, r = —1/3. 

The energies of the isolated levels can be obtained from 
equation cos(fcL) = 1, and the dispersed band energies 
can be found from equation cos(fcL) = [1 -I- 8cos(0)]/9. 
The ener gies of six lowest bands, see Fig. ||, are given by 
Eq. (|lOd|) and 



vrl : kL = arccos 



1 + 8 cos((/i) 



7r2 : kL = 27r — arccos 

7r3 : kL — 27r , 

7r4 : kL = 2n + arccos 

7r5 : kL = An ~ arccos 

7r6 : kL = An . 



l + 8cos(0) 
9 



l + 8cos(0) 

9 
1 -h 8 cos(0) 
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FIG. 3. The energies of calculated in the free electron ap- 
proximation PPV and PNV vr-bands. The only parameter in 
the calculations is perimeter of phenyl ring, L = 8.3A. 



The bands ttI - 7r4 are occupied, and bands 7r5 and 
7r6 are not. Really, we have 8 7r-electrons per unit cell 
and they occupy four lowest bands. Fortunately, our 
model gives no states near kL = 3tt/2, where the bands 
could be affected by the carbon chain periodic potential. 
However, the calculated energy gap between the top of 
the valence, 7r4, band, where kL — Sir — 2arcsin(l/3), 
and the bottom of the conductance, 7r5, band, where 
/cL = Stt — 2arcsin(l/3), is still smaller than the experi- 
mentally measured value 2.4-2.5eV. 

The symmetry of all obtained above states is essen- 
tially different from the benzene molecule states. It is 
because free electron wave vector k determines the co- 
ordinates dependence of the PPV wave functions and it 
changes continuously within the energy bands. The ex- 
plicit form of the calculated from Eqs. (aJq) PPV wave 
function amplitudes are given in Sec. |^ One can imme- 
diately see that direct optical transitions 7r4 — > 7r5 and 
7r3 — *■ 7r6 are dipole allowcdLl and their dipole moment in- 
terband matrix element is parallel to PPV chain. Direct 
optical transitions ttS —f 7r5 and 7r4 — > 7r6 are dipole for- 
bidden near the bands edges, (p = it. The dipole moment 
interband matrix element increases rapidly with transi- 
tion energy, and therefore, absorption bands correspond- 
ing to these processes also have to be experimentally ob- 
served. Our theory predicts following density of states 
peaks for four transitions: 1.4eV for 7r4 -^ ttS, 3.4 for 
ttS -^ 7r5, 4.4 for 7r4 -^ 7r6, and 6.4 for 7r3 -^ ttG. Three of 
these transtion energies are in qualitative agreement wi 
experimentally measured light absorption spectrurrjj 
which have features near 2.5eV, 3.7eV, 4.8eV, and 6.5eV. 

The similar picture is observed in PNV. It is mopj 
sophisticated conjugated system synthesized recently.EJ 
The levels of naphthalene molecule are also well described 
by the PFEO modelcl. In other words 7r-electrons move 
freely around naphthalene molecule, see Fig. || In this 
model naphthalene molecule is replaced by two bonds, 
one of them having length L/2 and another having length 
7L/6. Equations similar to Eqs. (y||§) can be written for 
this model of PNV and from the condition of their com- 
patibility we will obtain the following band structure: 

f 9 9 9 

cos(0) = <^ - cos{2kL) + - cos(5fcL/3) -I- - cos(4A:L/3) 

7 7 1-1 
-t- - cos(fcL) -I- - cos(2fcL/3) cos(fcL/3) \ 

8 8 8 J 

/ cos(fcL) -f cos(2fcL/3) -f 2 cos(fcL/3) -t- 1 . 

(12) 

This is algebraic equation of order five for e^^'^/^. Its 
solution leads to sequence of bands, which are shown 
on Fig. 0. PNV has 12 carbon atoms in unit sell and 
six lowest bands are occupied. There are no states near 
kL = 37r/2 in this system, however the calculated energy 
gap between valence, 7r6, and conductance, 7r7, bands is 
still smaller than measured valuai3 2.3eV. The conduc- 
tance band offset in possible PPV - PNV system is much 



smaller than the gap, that is in agreement with experi- 
ment. 



III. POLYMERS IN MAGNETIC FIELD. 

Let's consider experiment in which polymers are pre- 
pared in a two dimensional layer. All the phenyl rings 
of the polymer will be also lying in the same plane. If 
one applies magnetic field perpendicular to that plane, 
electron propagating around phenyl ring will acquire the 
additional phase, <&, which is equal to the magnetic flux 
via the ring divided by Tic/e. This additional phase can 
be easy introduced into Eqs. (||J|), and the determinant 
of the coefficients matrix becomes 



Q 



-, — ikL 



{r^t)H\{r + 2t)Hr-ty 



2r^ + 2rt - t^ + 2t^ cos($) 



4*2 cos(0) cos($/2) 



-ikL 



-2ikL 



(13) 



Again the equation Q = Q gives the dispersion relation 
for the electron states in the polymer. The states 7r3 
and 7r6 are not affected by the magnetic field, this can 
be explained by the properties of these states. In such 
a state electron is not only localized in thejphenyl ring, 
it is localized in one of the parts of the ringj. Therefore 
electron in such a state does not go around the phenyl 
ring and does not know about the magnetic field. 

The electron energy in other 7r-states is periodic func- 
tion of magnetic flux 



^4,5(0,*) = 



7r=F7((/),$) 



2mL2 



37rT7(0,*) 



2toL2 

1 4 
cos[7((/),$)] =---cos($) 

Q 

- -cos((/))cos($/2) . 
9 



(14) 
(15) 

(16) 



Here 7((/), $) > 0. The energy of 7r2 and 7r5 states de- 
creases in weak magnetic field (paramagnetic states) and 
the energy of ttI and 7r4 states increases (diamagnetic 
states). For example, the energy gap between the con- 
duction and the valence bands decreases with magnetic 
field, but the change is very small, since the area of the 
ring is very small, see Ref. |l^, Eq. 4-37. 

The above model does not take in to account the cur- 
vature of the polymers. If the polymer has the points 
of the self crossing the electron states can get additional 
magnetic moment that depends on the tunneling proba- 
bilities at such points. The magnetic moment of 7r-states 
without considering of this effect can be obtained from 



the derivatives of the energies with respect to magnetic 
field. 

The total magnetic moment of the polymer is given by 
the summation of the magnetic moments of all the occu- 
pied states. The diamagnetic susceptibility of the layered 
polymer material at zero temperature can obtained as a 
second derivative of the total energy of the system with 
respect to magnetic field. The approximate expression 
for the susceptibility per delocalized electron is 



A.= 



1 



P«/2 



X = N 






£2(0,$)+f4(0,$) 



(17) 



where S is phenyl ring area and N is number of 
monomers. In this expression we average susceptibil- 
ity over delocalized states in the energy interval corre- 
sponding to TT < fcL < Stt. One can check that the 
averaged susceptibility in interval < kL < tt is the 
same. It is convenient to express x in terms of obtained 
from PEEO theory benzene molecule susceptibility per 
electronlld, x^ — S^ e^ / {mc^ L^) . Numeric integration 
with respect to 4> gives x ~ O.TA^xo- Since PPV has six 
delocalized tt electrons per monomer, the PPV magnetic 
susceptibility per mole of monomers has to be equal 0.7 
of benzene molecule molar diamagnetic susceptibility: 

Xppv = 0.7xc6H6 (18) 



IV. SCATTERING OF FREE ELECTRON 
WAVES IN COPOLYMER. 

The last example for the application of our almost 
heuristic theory to the real physical problems is consid- 
eration of what will happen if we replace one monomer 
in a polymer by monomer of a different kind. Let's con- 
sider the copoly(p-phenylene vinylene - l,4-naplit|halene 
vinylene), which was used in recent experiments^. This 
system can be obtained formally from PPV if one re- 
places part of phenyl rings by naphthalene molecules. If 
the naphthalene concentration is small, we can consider 
each monomer of 1,4-naphthalene vinylene as an impu- 
rity or a structural defect in PPV chain. Therefore, such 
monomer partially reflects electrons propagating in PPV 
bands. 

We start the calculation of naphthalene vinylene 
monomer transmission coefficient with derivation of pure 
PPV wave functions. Three quantum numbers n, k, 
and s = sign(7) specify each PPV tt electron state, see 
Eq. (|lj) below. Therefore, the wave functions of all PPV 
TT states are given by the amplitudes A — Aa{n,s,K)^ 
B = BQ{n,s,K), C — Co(n, s, k), D — Z?o('t-, s, k), 
E — Eo{n, s, k), F = Fo{n, s, k), which can be calculated 
for given energy from Eqs. (||,||). Their explicit forms are 



Bn = 



Co — Eq — 



Do = Fo 



V2L|cosh(e)| 

1 

v/2L|cosh(e)|' 
i 



-«/2 



V4i|cosh(0| 



Jn/2+i/2 



jK/2-i/2 



v/4i|cosh(C)| 
cosh(^) = 3sin(7/2) , 
sinh(^) = 2\/2sin(K/2) , 

/o 

cos(7/2) = — — cos(k/2) 

O 

TT + K, > 0, K < 0, 

— TT + K, < 0, K > 0, 

koL = (2n + 1)tt + 7 , — tt < ^ < n 



(19a) 
(19b) 
(19c) 

(19d) 

(19e) 
(19f) 

(19g) 

(19h) 
(19i) 



where the wave function is normalized on one particle 
per unit cell of the polymer. In this system of equations 
free electron wave vector is also the function of quantum 
numbers k = ko(n, s, k). Vise versa, k determines the 
quantities k, 7, ^, n, and thus all of them are functions of 
energy. Dimensionless wave number k is more convenient 
than 0, which was used in Eqs. (11) and Fig. ^, because 
K goes to zero near the top of valence band and near the 
bottom of conductance band. There are two states for 
each pair of quantum numbers n and k, one state with 
positive 7, s = 1, and one state with negative 7, s = —1. 
For example the states in the conductance and valence 
bands are given by n = 1, with s = ±1 correspondingly. 
One can verify that the bands are parabolic near the 
edges, where |k| ^ 1. 

is not necessary, but introduced for 
When 7 is negative, the quantity ^ has 
imaginary part, which is equal to in. Since dimen- 
sionless wave number k have to be real Eq. 19g gives 
cos^(7/2) < 8/9. Therefore, the allowed values of 7 lie 
in the intervals — tt < 7 < —70 and 70 < 7 < tt, where 
7o = 2arcsin(l/3). The state is carrying electron, and 
flux is 



Equation |19g 
convenience. 



hk 
J ^ ^ tanh(^) , 
mL 



(20) 



which show the physical meaning of ^. 

Embedding of naphthalene vinylene monomer in to 
PPV chain result in scattering of electrons and holes 
with dimensionless wave number k. In order to calcu- 
late transmission, T, and reflection, R, coefficients we 
have to write equations similar to Eqs. (HJ^) for PNV 
free electron model, see Fig. ^ and substitute 

k = ko(n, s, k) , 

A = Ao{n, s, k) + i?Ao(n, s, —k) , 
A' = TBo{n, s, k) , 
B' = TBoin,s,K) . 



B = Bo{n, s, k) + RBQ{n, s, —k) 
B' = TB^{n, s, k) . 



(21) 



Here one have to use n = 1 and s = ±1 for conduction 
and valence band correspondingly. The further calcula- 
tions are straightforward, the qualitative results are the 
following. Near both conduction and valence bands edges 
transmission probability |Tp ss 0.8k^, and therefore it is 
proportional to the kinetic energy of quasiparticles. Far 
from the edges it becomes of the order of one. 



V. SUMMARY. 

We have developed a new free electron model of PPV, 
which allows to calculate band structure of polymer and 
hold wave functions of tt electrons in simple analytical 
form. We showed that the band structure contains only 
one parameter, which is the averaged distance between 
two neighboring carbons. The symmetry of PPV wave 
function allows direct optical transitions from ttS and 7r4 
states to 7r5 and 7r6 states, that explain four bands ob- 
served on absorption spectrum. The calculated energy 
gap between 7r4 and ttS states is not in agreement with 
experiment, but energies of three other transitions are 
in agreenie|nt. This suppjOijts. importance of many-body 
polaroni(£3 and excitonicE£lEZl effects for transitions near 
the energy gap. Similar models can be developed for 
other conjugated systems, for graphite planes or more 
complicated polymers. 

We have calculated exactly the magnetic field de- 
pendence of electron energies, which is due to pres- 
ence of phenyl rings. The diamagnetic susceptibility per 
monomer of PPV have been found approximately to be 
0.7 of the susceptibility anisotropy of isolated benzene 
molecule. We also considered scattering of electrons on 
polymer structural defects. We found that PPV elec- 
trons and holes transmission probability through 1,4- 
naphthalene vinylene monomer inserted in polymer is 
proportional to their kinetic energy. This is according 
to general property of scattering on potential occupiing 
final volume of space. We belive that proposed model 
have many potential applications. 
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